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ABSTRACT 

Context. We present time-resolved spectroscopy and photometry of HS 0218+3229, a new long-period cataclysmic vari- 
able discovered within the Hamburg Quasar Survey. It is one of the few systems that allow a dynamical measurement 
of the masses of the stellar components. 

Aims. We combine the analysis of time-resolved optical spectroscopy and _R-band photometry with the aim of measuring 
the mass of the white dwarf and the donor star and the orbital inclination. 

Methods. Cross-correlation of the spectra with K-type dwarf templates is used to derive the radial velocity curve of 
the donor star. An optimal subtraction of the broadened templates is performed to measure the rotational broadening 
and constrain the spectral type of the donor. Finally, an ellipsoidal model is fitted to the i?-band light curve to obtain 
constraints upon the orbital inclination of the binary system. 

Results. The orbital period of HS 0218+3229 is found to be 0.297229661 ± 0.000000001 d (7.13351186 ± 0.00000002 
h), and the amplitude of the donor's radial velocity curve is K2 = 162.4 ± 1.4 km s~^. Modelling the ellipsoidal light 
curves gives an orbital inclination in the range i — 59° ± 3°. A rotational broadening between 82.4 ± 1.2 km s~^ and 
89.4 ± 1.3 km s~^ is found when assuming zero and continuum limb darkening, respectively. The secondary star has 
most likely a spectral type K5 and contributes ~ 80 — 85% to the _R-band light. Our analysis yields a mass ratio of 
0.52 < q < 0.65, a white dwarf mass of 0.44 < Mi{Mq) < 0.65, and a donor star mass of 0.23 < M2{Mq) < 0.44. 
Conclusions. We find that the donor star in HS 0218+3229 is significantly undermassive for its spectral type. It is 
therefore very likely that it has undergone nuclear evolution prior to the onset of mass transfer. 

Key words, accretion, accretion discs - binaries: close - stars; individual: HS 0218+3229 ~ novae, cataclysmic variables 



1. Introduction 

Testing the current evolutionary models of the present-day 
Galactic population of cataclysmic variables (CVs) is not 
possible without a significant sample of the actual masses 
of the white dwarfs (i.e. the accretors) and the late-type 
companion stars (i.e. the donors) that form this sort of 
binary system. But accurate measurements of the compo- 
nent masses in CVs are, however, not an easy task. The 
donor stars are usually veiled by the extreme brightness 
of the accretion structures, especially the accretion disc. 
Furthermore, the white dwarfs are also very often hid- 
den from view in CVs with orbital periods above the 2- 
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3 hour period gap, so optical studies aimed at dynamical 
mass measurements are scarce and usually make dubious 
assumptions. A common practice has been to assume that 
the amplitude of the radial velocity curve of the wings of 
the disc emission lines can be regarded as a reasonable es- 
timate of the radial velocity amplitude of the white dwarf 
(i^i), but this assumption is usually far from being true 
(see e.g. Shaffer et al. 1995). Therefore, time-resolved stud- 
ies of CVs in which the donor star and/or the white dwarf 
are exposed are chiefly important to our understanding of 
CV evolution. 

HS 0218+3229 was discovered during our large-scale 
search for CVs in the Hamburg Quasar Survey (HQS), 
based on their spectroscopic properties or, more specifi- 
cally, the presence of strong emission and/or absorption 
lines (|Gansicke et all l2002b[ i. Fifty-three new CVs were 
found within this project, and the orbital period has been 
determined for the majority of them. The main result 
is that most of the HQS CVs have been found to have 
orbital periods above the period gap. Those include rarely 
outbursting dwarf novae, such as GYCnc (HS 0907+1902, 
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Gansicke et aII[20Qg) or RX J0944.5+0357 (HS 0941+0411, 



Mennickent et all l2002f ): magnetic CVs with relatively 
weak X-ray emission s uch as IRXS J062518 . 2+ 73343 3 

[Araujo-Betancor et al.1 



emission 
(HS 0618+7336, 

RXJ1554.2+272 1 (HS 15 52+27 30, 



, 12001 

iThorstensen fc FentonI 
and HS 0943+1404 
a n umb er of SWSextantis 
120011 : iRodriguez-Gil et all 
12004 I2007fl : on e of the youngest pre -CVs known 
CHS 1 857+5144. lAungweroiwit et all 120071): long-period 
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stars 



(ISzkodv et al 



■ .20041) . 
2005h : 



CVs dAungweroiwit et al.1 2005l i: and several dwarf novae 



(jAungweroiwit et al 200' 



This paper presents the dynamical determination of the 
stellar masses in the cataclysmic variable HS 0218+3229. 
We present the data in Sect. [21 and an identification spec- 
trum in Sect. [3] The orbital period is measured from the 
optical photometry in Sect. 3] before analysing the time- 
resolved spectroscopic data in Sect. [5] Modelling of the el- 
lipsoidal modulation is carried out in Sect. [6] Finally, the 
system parameters of HS 0218+3229 and the overall discus- 
sion are presented in Sect. [71 



HS0218+3229 .. ^j,':;. 7, , - M-+_ 

Fig. 1. 10' X 10' finding chart of HS0218 obtained from 
the Digitized Sky Survey 2. The coordinates of the new 
CV are a(J2000) = 2h21'"33.50", (5(J2000) = +32°43'23.8". 
The star 'CI' has been used as comparison for the i?-band 
differential photometry. 



2. Observational data 

2.1. Photometry 

Astrophysikalisches Institut Potsdam. Time series R- 
band photometry of HS 0218+3229 (hereafter HS0218) 
was obtained on 3 nights during the period 2001 
January-February using the 0.70 m telescope of the 
Astrophysikalisches Institut Potsdam (AIP), Germany. 
The images were obtained with the camera equipped 
with a SITe 1024x1024 pixel CCD detector and reduced 
in a standard way with MIDAS. Point spread function 
(PSF) photometry was done with a Perl Data Language 
pipel ine based on the DoPhot package (M ateo fc Schechteij 
Il989f) . The differential magnitudes of HS 0218 were derived 
relative to the comparison star labelled 'CI' in Fig. [2 
and were converted into apparent i?-band magnitudes 
using the USNO-A2.0 magnitude of the comparison star, 
Rci = 13.74. The main source of uncertainty in this 
conversion is the uncertainty in the USNO magnitudes, 
which is typically ~ 0.2 mag. 

Fred Lawrence Whipple Observatory. Further time-resolved 
i?-band photometry of HS0218 was secured on 2005 
October 8-11 using the 1.2 m telescope at the Fred 
Lawrence Whipple Observatory (FLWO) in Arizona. The 
images were obtained with the KEPLERCAM mosaic cam- 
era which consists of an array of four 2048x2048 pixel CCD 
detectors. Only a small window on CCD #4 was read out in 
order to minimise the dead time between images. Aperture 
photometry was performed on the reduced images using 
IRAlQ- An optimal aperture radius of 1.5 times the full- 
width at half -maximu m (FWHM) of the typical seeing pro- 
file was used ()Navloill99 8). The light curves were computed 
in the same manner as described for the Potsdam data. A 
brief summary of the observations is given in Tabled 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under contract 
with the National Science Foundation. 



Table 1. Log of observations 



UT Date Coverage 


Filter/Grating 


Exp. 


# Frames 


(h) 




(s) 




Calar Alto 2.2 m, CAFOS spectroscopy 




2000 Sep 20 


B-200/R-200 


600 


1/1 


AIP photometry 








2001 Jan 08 4.09 


R 


60 


223 


2001 Jan 16 6.11 


R 


60 


330 


2001 Feb 15 4.78 


R 


60 


251 


Calar Alto 3.5 m, TWIN spectroscopy 




2002 Oct 28 6.63 


T08/T01 


600 


38 


2002 Oct 29 3.10 


T05/T06 


600 


18 


FLWO 1.2 m photometry 






2005 Oct 08 5.65 


R 


60 


244 


2005 Oct 09 7.52 


R 


60 


377 


2005 Oct 10 7.17 


R 


60 


372 


2005 Oct 11 7.97 


R 


60 


413 



2.2. Optical spectroscopy 

During an identification run of HQS CV candidates at 
the 2.2 m Calar Alto telescope on 2000 September 20, we 
obtained a pair of blue/red spectra of HS0218 with the 
CAFOS spectrograph (Table [T]). We used the B-200 and 
R-200 gratings in conjunction with a 2" slit, which pro- 
vided a spectral resolution of ~ 10 A (FWHM). A standard 
reduction was performed using the CAFOS MIDAS quicklook 
package. The detection of strong Balmer emission lines and 
absorption lines of a late-type star confirmed the CV na- 
ture of HS 0218 and encouraged us to make the follow-up 
observations described in what follows. 

Time-resolved spectroscopy of HS0218 was obtained 
with the double- armed TWIN spectrograph on the 3.5 m 
telescope in Calar Alto on 2002 October 28-29 (Table [J). 
On the first night, a total of 38 blue and red spectra 
were acquired using the T08 and TOl gratings, respec- 
tively, and a 1.2" slit. The wavelength ranges AA3800 — 5560 
and AA5800 — 7390 were sampled at 2.2 and 1.6 A reso- 
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Fig. 2. Discovery spectrum of HS0218 obtained with the 
2.2 m telescope at Calar Alto in September 2000 (solid line). 
The spectrum of a K4 V star (HD 5351) is shown as a 
dotted line, shifted down by 0.5 units for clarity. The tem- 
pl ate spectrum was ex tracted from the spectral catalogue 
of iJacobv et al.l (|1984D . The absorption-line spectrum of 
HS0218 is characteristic of a mid-K star. 

lution (FWHM), respectively. The gratings T05 and T06 
and the same slit width were selected for the second night. 
This granted access to the ranges AA3985 — 5050 and 
AA5910 - 6990 at 1.1 and 1.2 A resolution (FWHM) in the 
blue and the red, respectively. Spectra of the G191-B2B 
flux standard were taken to derive the instrument response 
function. Both the identification spectra and the spectra 
obtained on the first night showed a red continuum and ab- 
sorption features originating in the photosphere of a mid-K 
secondary star. Hence, we also took spectra of six K stars 
of different spectral sub- type during the second night. For 
wavelength calibration we acquired spectra of a HcAr lamp 
throughout each night. 

Prior to optimal extraction of the individual spectra 
(|Hornelll986n .the raw images were bias-subtracted and flat- 
field corrected. Sky emission was removed during the opti- 
mal extraction of the spectra. The pixel-wavelength data 
were well fitted by a fourth-order polynomial, with an rms 
always smaller than one tenth of the spectral dispersion in 
all cases. The wavelength scale for each spectrum was de- 
rived by interpolating between the two nearest arc spectra. 
The accuracy of our wavelength calibration was improved 
by using the night-sky lines to correct for any zero-point 
offset. 

The reduction steps and one-dimensional spectrum ex- 
traction were performed within IRAF, whilst the wavelength 
calibration was performed in MOLL\Q. 

3. Identification spectroscopy 

One of the remarkable features of the identification spec- 
trum of HS0218 (Fig. [2]) is its red continuum, which re- 
veals the dominance of the companion star. The absence 
of strong TiO absorption bands suggests a spectral type 

^ Written by T. R. Marsh, 
http : / /www. Warwick, ac.uk /go /trmarsh 



To 


Cycle 


O-C 


(HJD - 2450000) 




(s) 




1918.39911 


-5836 


239 


1926.27401 


-5809.5 


94 


1926.42749 


-5809 


514 


1956.30491 


-5708.5 


1018 


2576.44547 


-3622 


-1496 


2577.33716 


-3619 


-1495 


3652.88176 


-0.5 


156 


3653.77448 


2.5 


245 


3653.91820 


3 


-177 


3654.67436 


5.5 


953 


3654.80932 


6 


-226 


3654.95944 


6.5 


-96 


3655.70444 


9 


70 


3655.85465 


9.5 


208 



earlier than M. In fact, the broad absorption observed at 
- 5200 A (produced by MgH at 5180 A, a TiO band at 
AA4954 — 5200, and a jump in the continuum due to the 
Mg I triplet at AA 5168 — 5185) constitutes a clear signature 
of a K-type star (jMcClintock fc Remillardi ri990. and ref- 
erences therein). In Sect. 15.21 we present a more accurate 
determination of the secondary star spectral type based on 
measurements of the rotational broadening of its absorption 
lines. 

4. The orbital period of HS 0218 

The i?-band fight curves of HS 0218 (Fig.g]) exhibit a clear 
quasi-sinusoidal modulation with variable minima. Deeper 
and shallower minima alternate in a periodic basis every 
^ 3.5 h. This behaviour, together with an optical spectrum 
dominated by the donor star (Fig. [2]) , make us identify this 
variation with a classical ellipsoidal modulation resulting 
from the changing projected area of the Roche lobe with 
orbital phase. The deeper minima in the ellipsoidal light 
curve are caused by the stronger gravity darkening on the 
hemisphere facing the white dwarf, so they must occur at 
orbital phase 0.5. Therefore, the actual orbital period of 
HS 0218 must be twice the separation between consecutive 
minima, which is confirmed by the radial velocity curve of 
the donor star (Sect. [5T|) . 

As the photometry spans much longer than the spec- 
troscopy, the analysis of the light curves provided a more 
accurate orbita l period determination. Th e analysis-of- 
variance (AOV, ISchwarzenberg-Czernvl[T989t ) periodogram 
computed from our light curves is presented in Fig. [H 
The highest peak is centred at i/ ~ 6.73 d~^ (~ 3.57h), 
whilst the second highest lies exactly at half that frequency 
(ri ~ 3.36 d^^; Porb — 7.13h), which is the actual orbital 
frequency. 

In order to improve the accuracy of the orbital pe- 
riod we calculated a linear ephemeris from the times of all 
available photometric minima derived from Gaussian fits. 
Additionally, two measurements of the instant of zero phase 
(i.e. inferior cojunction of the secondary star, Tq) were ob- 
tained from sine fits to the secondary star's radial velocity 
curves presented in Sect. 15. li which were also included in 
the calculations. The measurements of the times of phase 
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Fig. 3. i?-band light curves of HS 0218+3229 obtained at the AIP [left) and the FLWO [right). 



(and also 0.5 for the ellipsoidal light curves) are presented 
in Table [2l The resulting orbital ephemeris is: 

To(HJD) = 2 453 653.028599(3) + 0.297229661(1) x S , (1) 

where the numbers in parentheses quote the uncertainty in 
the last digit. 

Figure [4] also shows all the i?-band data folded on the 
orbital period using the ephemeris given in Eq. [1] The el- 
lipsoidal modulation is apparent. The shape of the folded 
light curves, however, shows significant deviations between 
the AIP and FLWO data. The 2001 January 8 light curve 
clearly displays the phase-0.5 minimum and phase-0.75 
maximum shifted by ~ —0.05 cycle with respect to the 
FLWO light curve. This cannot be the eff'ect of an inac- 
curate orbital period determination as the phase-0 minima 
and phase 0.75 — branches match very well in both sets of 
data. In addition, the 2001 January 16 AIP light curve ex- 
hibits a broader phase-0 minimum, with the phase — 0.25 
rising branch lagging that of the FLWO and almost equal 
primary and secondary minima. This behaviour indicates 
that the ellipsoidal modulation is actually contaminated by 
other sources of variation, likely originating in the accretion 
disc and/or the donor star itself (e.g. irradiation). 

On the other hand, the i?-band light curves of 
HS0218 also display higher frequency activity (see e.g. the 
2005 October 10 light curve in Fig. [3]). To analyse these 
variations we first subtracted the ellipsoidal modulation (by 
means of a smoothed, Fourier series-fitted version of the 
curves) and then computed an AOV periodogram from all 
the detrended light curves. The combined power spectrum 
does not show any predominant frequency. However, the pe- 
riodogram computed from the 2005 October 9 data alone 
(not shown) exhibits a prominent peak at ~ 37.4 d~^, 
which corresponds to a period of P ~ 38.9 min (= 0.027 d), 
approximately 10% of the orbital period. The periodograms 
calculated from the other individual light curves (with the 
exception of 2005 October 8) exhibit their strongest peaks 



at close but lower frequencies. This might suggest the pres- 
ence of quasi-periodic oscillations in HS0218, but a much 
more longer light curve coverage would be needed to draw 
any firm conclusions. 

5. Time-resolved spectroscopy 

5.1. The radial velocity amplitude of the donor star 

The radial velocity curve of the secondary star in HS0218 
was measured from the red spectra using the method of 
cross-correlation with the spectrum of a template star 
(iTonrv fc David [19791. Five template stars with spectral 
types ranging from KO V to K5 V were acquired on the night 
of 2002 October 29 to extract the radial velocities. Prior 
to the cross-correlation, the spectra were re-sampled onto 
a common logarithmic wavelength scale and normalised by 
means of a low-order spline fit to the continuum. Individual 
velocities were then extracted by cross-correlation with 
each template star in the range AA6090 — 6520. We chose 
this spectral interval because it contains lines useful not 
only for radial velocity measurements, but also for spectral- 
type classification and rotational broadening measurements 
(see Sect.lOl)- 

A first look at the velocities indicates a period of ~ 0.3 d, 
which confirms the ellipsoidal nature of the i?-band light 
curves (Sect.S]). We therefore performed least-squares sine 
fits to the radial velocity data with the period fixed at the 
photometric orbital period. The resulting fit parameters are 
summarised in Table [3) Although all the values are very 
similar, we adopt the K5 V parameters of the radial velocity 
curve since the spectral type of the donor star is shown to 
be most likely K5 (see Sect. 15.2]) : 

K2 = 162.4 + 1.4 km s-i 
72 = -54.2 + 2.9 km s^^ 

All quoted uncertainties are 1-a. The phase-folded radial 
velocity curve is shown in Fig. [5l 



p. Rodriguez-Gil et al.: An evolved donor star in HS 02184-3229 
Table 3. Radial velocity parameters (errors at la). 



5 



Template 


Spectral 


7 


72 '"^^ 




To 






type 


(km s"^) 


(km s~^) 


(km s-i) 


(HJD - 2452576) 


(d.o.f=53) 




HD 10780 


KOV 


-56.2 ± 1.1 


-53.5 ± 1.2 


161.6 ± 1.5 


0.4457 ± 5x10"" 


9.1 


HR 10476 


KIV 


-20.9 ± 1.1 


-54.7 ± 1.1 


161.7 ± 1.5 


0.4455 ± 5x10^" 


9.4 


HR 16160 


K3V 


-76.4 ± 1.0 


-50.7 ± 1.2 


162.3 ± 1.4 


0.4457 ± 4x10"* 


10.1 


HD 6660 


K4V 


-59.7 ± 1.0 


-54.9 ± 1.6 


162.2 ± 1.4 


0.4456 ± 4x10"" 


13.5 


HD 10436 


K5V 


-2.8 ± 1.0 


-54.2 ± 2.9 


162.4 ± 1.4 


0.4457 ± 4x10"'' 


11.9 



Systemic velocity after adding the radial velocity of the template star as found in Gontcharov (2006). 



Table 4. Spectral classification and rotational broadening. 



Template 


Spectral 


vsini (km s"^) 


./ 


xl 


V sin i '''' (km s """) 




type 


(m = 0.5) 




(d.o.f.=496) 


[/i = 0.0, /i at continuum] 


HD10780 


KO V 


89.5 ± 1.2 


1.20 ± 0.01 


2.1 


[84.7 ± 1.8, 91.4 ± 2.0] 


HD 10476 


Kl V 


89.6 ± 1.2 


1.17 ± 0.01 


1.8 


[84.8 ± 1.6, 91.2 ± 4.3] 


HD16160 


K3 V 


89.0 ± 1.2 


0.97 ± 0.01 


1.0 


[84.3 ± 1.2, 91.1 ± 1.3] 


HD 6660 


K4 V 


87.7 ± 1.2 


0.728 ± 0.007 


1.0 


[83.2 ± 1.1, 90.0 ± 1.1] 


HD 10436 


K5 V 


86.9 ± 0.6 


0.785 ± 0.008 


0.9 


[82.4 ± 1.2, 89.4 ± 1.3] 



Uncertainties are 1-cr corresponding to Xmin + 1 (Lampton, Margon & Bowyer 1976). 

Uncertainties are 1-a and were estimated using the Monte-Carlo method by simulating a total of 1000000 
copies of our target spectrum using the bootstrap technique outlined in Steeghs fc Jonker (2007). 



5.2. Spectral type and rotational broadening 

To determine the spectral type of the donor star and the 
rotational broadening of its absorp t ion li nes, we used the 
technique outlined in iMarsh et al.l (|1994| ). It is based on 
the search for the lowest residual obtained when subtract- 
ing a set of templates from the Doppler-corrected, average 
spectrum of the target. 

We proceeded as follows: we used the spectra acquired 
during the second night, as both target and template spec- 
tra were acquired with the same instrumental setup. This is 
essential to derive a reliable value for the rotational broad- 
ening (v sin i) , in particular when v sin i is of the same order 
as the instrumental resolution. First, the 18 target spectra 
were Doppler-corrected to the rest frame of the donor star 
by using the radial velocity parameters derived in Sect. 15. ll 
Next, we produced an average spectrum after assigning dif- 
ferent weights to the individual spectra to maximise the 
signal-to-noise ratio of the sum. The template spectra were 
broadened from 70 to 110 km in steps of 1.0 km s~^ 
through convolution with the spherical rotational profile of 
Gray (1992) with a linearised limb-darkening coefRcient of 
0.5. Each broadened version of a template spectrum was 
scaled by a factor / (representing the fractional contribu- 
tion of light from the secondary star) and subtracted from 
the target Doppler-corrected average. Then a test on the 
residuals was performed in the range AA6090 — 6520 and the 
optimal values of / and v sin i were provided by minimising 

Besides the instrumental broadening, the secondary 
star's photospheric lines are also smeared due to the or- 
bital motion of the donor star during the length of the 
exposure. With 600-s exposures, the orbital smearing is 



< texp^K2 = 24 km s"^. This effect can be corrected pro- 
ducing a spectral template by averaging 18 smeared copies 
of a template spectrum with the same weights as used to 
average the spectra of the target. We performed this cor- 
rection, finding no significant difference in the values cal- 
culated for the rotational broadening. 

The minimisation of xt (Table 2]) shows that the spec- 
tral type of the secondary star in HS 0218 is later than Kl. 
Also, a KO-1 V donor star is rejected as the templates give 
/ > 1. Similar results for the spectral classification were 
obtained when repeating the above procedure on the 38 
target spectra acquired during the first night. Additional 
constraints on the spectral type of the donor star are pro- 
vided by the infrared colours of HS0218 after correcting 
for reddening and assuming the same fdisc at each infrared 
bandpass. Using 2MASS magnitudes we derive {J — H, iJ — 
K, J -K)n = (0.63 0.09 0.72) based on E{B-V) = 0.071 
mag (jSchlegel et al.|[T998l) . The uncertainties in the infrared 
colours are 0.04 mag. The above colours reject an M-type 
dwarf {J -K > 0. 86, H - K > 0.165) a nd support a K4- 
K5 dwarf (see e.g. iBessell fc Br"ettlll988f ). For comparison, 
a K5 dwarf has colours of (0.61, 0.11, 0.72). Further visual 
inspection of the average spectra and templates supports 
a K5 spectral type based on the strength of the weak TiO 
bands in the spectrum and by comparing the relative inten- 
sities of the absorption lines in the interval AA6400 — 6530. 
Hereafter, in accordance with the above analysis, we adopt 
a K5 spectral type for the donor star. 

A rotational broadening measurement of 86.9 ± 0.6 km 
s"^ was obtained from the K5V template. However, the 
value of the limb-darkening coefficient for the absorption 
lines is unknown and this limits our ability in evaluating 
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0.1 



Orbital phase 

Fig. 4. Top main panel: analysis-of-variance (AOV) pe- 
riodogram computed from the i?-band photometric data 
(Fig. [3]). The periodogram shown on top was constructed 
from a fake sine wave with a frequency of 20 sampled at 
the actual data points. The structure of both periodograms 
are nearly the same. Bottom panel: the i?-band photometric 
data folded according to the ephemeris given in Eq. [1] The 
solid line is the best ellipsoidal fit to the curve (see Sect. [6]). 
No phase binning has been applied. A whole orbit has been 
duplicated for clarity. 



V sin I. Limb-darkening coefficients for the absoption lines in 
late-type stars are expected to be smaller than those for the 
continuum (Collins & Truax 1995). Therefore, we repeated 
the above analysis for both a null linearised limb-darkening 
coefficient and for the value expected for th e continuum. 
The la tter was selected for each template from lClaret et al.l 
()1995( ). We provide the results for each template in Table 
3 (column 6). As one can see, our uncertainty in the limb- 
darkening leads to values of vsini of 82.4 — 89.4 km s~^ 
(K5V template). 



5.3. The radial velocity amplitude of the white dwarf 

We measured the radial velocity curve of the Ha line 
wings by using the double-Gaussian technique of 
[Schneider fc Yound (|1980D . In an ideal world, with an ax- 
isymmetric accretion disc dominating line emission, the 
line wings would provide a good estimate of Ki. We ob- 
tained radial velocity curves of the Ha wings for a num- 
ber of Gaussian (FWHM = 200km s~^) separations be- 
tween 600 and 2900km s^^ in 100km steps. The curve 
for a 1300km s^^ separation gave the best results and is 
shown in Fig. \5\ The wings are clearly delayed by ~ 0.1 
cycle, but their 7- velocity of —57.6 ± 0.8 km s^^ (derived 
from a sine fit) is very similar to that of the donor star 
(72 = —54.2 ± 1.4km s~^). The sine fit gives an amplitude 
of the wings radial velocity curve of 100.3 ± 1.2km s~^, 
which provides a mass ratio oi q = K1/K2 = 0.618 ± 0.009 




Orbital phase 

Fig. 5. Radial velocity curves of the secondary star (solid 
dots) and the Ha line wings (open circles) in HS0218. The 
radial velocities of the donor star were extracted by cross- 
correlation with a K5 V template star. Solid curves are the 
best sine fits to the radial velocities. The horizontal dashed 



line marks 72 — 
plotted twice. 



54.2 km s . The orbital cycle has been 



if assumed to be the radial velocity amplitude of the white 
dwarf. This mass ratio is equal to the measured from the 
rotational broadening of the donor star absorption lines. 
Despite the phase shift that the Ha radial velocities show 
(Fig. ini) , the measured Ki value seems to be reliable in this 



6. Ellipsoidal light curve modelling 

To model the ellipsoidal variation we used a binary code 
which uses full Roche geometry and a black body model 
to account for the emission of the distorted companion. 
In order to carry out the numerical integration the donor 
star surface is divided in small triangular tiles of equal area 
which cover completely the Roche-lobe filling star. In this 
model, the disc contribution to the total light (/disc = 1 — /) 
is independent of the orbital phase and is scaled to a cer- 
tain factor of the mean flux from the companion. The limb 
darkening coefficients are obtained from Al-Naim iy ( 19781) 
and t he /3 gravity da rkening expo nent is fixed at 0.08 (jLucvl 
|1967[ ). According to iGravl (|1992D the effective temperature 
of a K5V star is 4557 K, which was adopted as the polar 
temperature. The synthetic ellipsoidal modulation should 
only depend on the orbital inclination angle, the mass ratio 
{q = M2/M1, where M2 and Mi are the masses of the donor 
star and the compact object, respectively), and the relative 
disc contribution to the light curve (1 — /). In Sect 15.21 we 
obtained 1 — / ^ 0.22 by measuring the absorption lines in 
the range AA6090 — 6520. This spectral region is contained 
in the R band and we can therefore use this value in our 
fits. However, we think that the contamination in the shape 
of the light curve may be due to night-to-night variations 
of the disc contribution (see Sect. S]). Hence, the value ob- 
tained during the spectroscopic campaign may not be the 
correct one for the photometric data. 

In order to establish a conservative lower limit to the 
inclination angle (i) of the system we have fitted the light 
curve by using 1 — / = 0, obtaining i > 50°. Since the light 
curves do not show eclipses we find an upper limit of j < 



p. Rodriguez-Gil et al.: An evolved donor star in HS 0218+3229 



7 



62° from elemental trigonometry involving the maximum 
disc radius allowed by tidal interaction with the donor star 
(jPaczvh ski 1977) and the Roche lobe radius of the latter. 
Thus, we obtain 50° < j < 62° independently of /disc- The 
best fit to the ellipsoidal light curve was finally found for 
/disc — 0.15, q = 0.62, and i ~ 59 ±3°. We have overplotted 
this synthetic model to the actual folded ellipsoidal curve 
in Fig. m 



7. Discussion 

7.1. Parameter estimates 

Assuming that the companion star is synchronised with 
the binary motion and fills its Roche lobe, it is possi- 
ble to calculate the mass ratio through the expression 
vsmi = 0.462 K2 q^/^{l + q)^^^ (e.g. Wade & Home 1988). 
Thus we derive 0.52 < g < 0.65 from K2 = 162.4 ± 1.4 km 
and the 99% confidence range for the rotational broad- 
ening (82.4 - 2.33a km s"^ < vsini < 89.4 -I- 2.33cr km 
s""'^). This mass-ratio range implies a primary velocity am- 
plitude of 84 km s^i < Ki^ qK2 < 106 km s^^. The value 
of 100.3 ± 1.2 km s"^ measured from the Ha line wings lies 
almost in the middle of this interval. 

The velocity amplitude and orbital period imply a mass 
function of /(Mi) = PK^/2ttG = 0.132 ± 0.001 M© = 
Misin'^i/(1 -I- q)^. Combining the values for /(Mi) and 
our extreme limits on q and i, we constrain the masses of 
the compact object and the companion star to be in the 
range: 0.44 < Mi < 0.65 Mq and 0.23 < M2 < 0.44 Mq. 
Table [5] gives the determined system parameters for the 
extreme values of the inclination given by the ellipsoidal 
modelling. The uncertainties in the masses for a given in- 
clination are dominated by the uncertainty in the value of 
usini. For comparison, using the resulting value obtained 
for a limb-darkening coefficient of 0.5 {vsini = 86.9 ± 0.6 
km s-i), we derive: q = 0.61 ±0.01, Mi = 0.54 ±0.03 Mq, 
M2 = 0.33 ± 0.02 M© and Xi = 99 ± 1 km s'^ To deter- 
mine the above 1-a uncertainties a Monte Carlo approach 
was used in which we draw one million random values for 
the observed quantities. We treated K2 and vsini as be- 
ing normally distributed about their measured values with 
standard deviations equal to the uncertainties on the mea- 
surements. For i the distribution was taken to be uniform. 
Given the stellar masses and orbital period, we can con- 
strain the semi-major axis (a) of the orbit using Kepler's 
Third Law. the radius of the secondary star {R2) using 
Egglcton's (1983) expression for the effective radius of the 
secondary Roche Lobe, and the luminosity of the secondary 
star (L2) using Stefan-Boltzmann's Law and the 4557 K ef- 
fectiv e temperature inferred from the spectral type (|Gravl 
|1992|) . These quantities are listed in Table [H] 

To constrain the distance to HS0218 we have adopted 
the luminosity of the donor star given in Table [51 which, 
combined with the bolometric correction and the {V — R) 
colour of a K5 V star, yield an absolute magnitude of Ma — 
5.96 — 5.65. We have calculated a dereddened magnitude of 
the secondary star of R — 15.65 — 15.72 using the apparent 
magnitude of HS 0218 (i? = 15.6), Ar = 0.188 and a - 15- 
20% veiling from the accretion disc (Sect.O. Therefore, the 
distance modulus relation provides a distance to HS0218 
of 0.87 — 1.0 kpc. Note that this distance estimate must be 
considered as approximate given the uncertainties in our 



Table 5. Summary of derived system parameters. 



Parameter 



56° 



: 62° 



Ml (Mq) 
M2 (Mq) 

R2 (Rq) 
a (Rq) 

L2 (Lq) 



0.54 - 0.65 0.44 - 0.54 

0.28 - 0.44 0.23 - 0.36 

0.57 - 0.67 0.53 - 0.63 

1.76 - 1.93 1.65 - 1.81 

0.12-0.17 0.11-0.15 



photometric calibration based on USNO A2.0 magnitudes 
of field stars. 



7.2. An evolved donor star 

As we have demonstrated, the companion star in HS0218 
is most likely of spectral type K5. A main-seq uence K5 
star has a nominal mass of 0.68 Mq (lGravlll992D . The dis- 
agreement with the mass constraints obtained above is eas- 
ily explained if the secondary star is undermassive for its 
spectral type, that is, has undergone nuclear evolution be- 
fore the onset of mass transfer. In fact, with a mass of 
~ 0.33 Mq its nuclear timescale is much greater than the 
Hubble time. Previous wor k on donor stars in long -period 
CVs (see e.g. figure 5 in iBeuermann et al.l Il998[ ) shows 
that most of them at long orbital periods are underlumi- 
nous, and likely evolved. In Fig. [5] we plot evolutionary 
sequences for an initial donor star mass of 1.0 Mq and 
mass transfer onset at several central hydrogen abundances 
(jBaraffe fc Kolbll2000D . Note that this is only for illustrative 
purposes and the actual initial donor mass in HS0218 can 
be different. The mass of the donor star in HS 0218 is much 
lo wer than predicted by an y of the sequences considered. 

ISchenker fc Kinj (|2002l ) show the evolution of the donor 
mean density as a function of its mass (their figure 2) for 
both the ZAMS and evolved cases. With an orbital pe- 
riod of 7.13 h, HS0218 has a mean density of 2.1 g cm~^, 
which results in a donor mass of ^ 0.3 Mq for the case of 
strong thermal-timescale mass transfer. Moreover, accord- 
ing to Schenkex_et_al, ( 2002) (their figure 6), a donor star 
of spectral type K5-K7 is expected for a mass of ^ 0.3 Mq. 
Both predictions for the case of an evolved donor agree well 
with our measurements. 

The binary system, therefore, formed with a compan- 
ion star more massive than the white dwarf. In this sce- 
nario {q — M2/M1 > 1), and as a consequence of mass 
transfer towards the white dwarf, the Roche lobe of the 
donor star (initially the most massive component) shrinks 
faster than the thermal equilibrium radius. The secondary 
star thus tends to expand beyond the Roche lobe to reach 
the thermal equilibrium radius, which triggers an episode 
of thermal-timescal e mass transfer until the m ass ratio 
reverses (see e.g. King_et_al.| ■200l ISchenker et' al. 200j). 
Beyond this point, normal mass transfer starts due to ei- 
ther angular momentum loss (i.e. secondary Roche lobe 
shrinkage) of nuclear evolution of the donor star (i.e. in- 
crease of the thermal equilibrium radius) and the system 
becomes a normal CV. The fact that the measured white 
dwarf ma ss (Mi ^ 0.54 Mq) is lo wer than the average value 
for CVs (jSmith fc Dhillonlll998l) indicates that the trans- 
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5 E) 7 

OrbiLal period (h) 

Fig. 6. Spectral type, mass, and radius as a function of or- 
bital period for four different evolutionary sequences with 
different initial central hydrogen abundance (^c) at the on- 
set of mass transfer. The initial donor mass is M2 = 1.0 M0 
and the mass transfer rate is M — 1.5 x 10~^ M0 yr^^. 
Solid line: — 0.7 (initially unevolved donor); dotted 
line: Xc = 0.55; dashed line: Xc — 0.16; long-dashed 
Hne: Xc = 0.05. The location of HS0218 is indicated by 
the large black dot. E volutionary sequences adapted from 
iBaraffe fc Kolbl (|2000h . 



ferred donor material during the thermal-timescale stage 
has been blown off the binary system. 

Nuclear evolution of the donor star prior to normal CV 
life should reflect in the chemical composition of the ac- 
creted material. Because of CNO processing of the ma- 
terial, carbon must be significantly depleted and nitrogen 
subsequently enriched. Unusually low Civ/Nv ultraviolet 
line flux ratios have been observed in a number of CVs 
(jGansicke et al.l [20031 ). suggesting that a significant frac- 
tion of the present-day CV population is descendant from 
systems that have passed through a thermal-timescale mass 
transfer phase. HS0218 is therefore a promising target for 
ultraviolet spectroscopy to check whether N enrichment has 
taken place. 
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